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ON SOME EXTREMALITIES IN THE APPROXIMATE INTEGRATION 



SZYMON WASOWICZ 



ABSTRACT. Some extremalities for quadrature operators are proved for convex functions 
of higher order. Such results are known in the numerical analysis, however they are of- 
ten proved under suitable differentiability assumptions. In our considerations we do not 
use any other assumptions apart from higher order convexity itself. The obtained inequal- 
ities refine the inequalities of Hadamard type. They are applied to give error bounds of 
quadrature operators under the assumptions weaker from the commonly used. 



1. Introduction 
In the theory of convex functions the Hermite-Hadamard inequality 

which holds for convex functions (and, in fact, characterizes them), plays a very important 
role. The first inequality follows by the existence of a support line for / at the midpoint, 
while the second one can be obtained using the fundamental property of convexity stat- 
ing that a graph of a convex function / lies on [a, b] below the chord joining the points 
(a, /(a)), (6, /(&)). We have also the following 

Observation 1. If a real function f is convex on an interval [a, b] then 

for any N 6 N, £i, . . . , £jy 6 [a, b] and Ai, . . . , A^v ^ with Y^iLi A< = 1 such that 
S^ N \ C — 2d± 

Proof. The first inequality is an immediate consequence of convexity, the second one we 
prove similarly to the second inequality of (Q~|). Q 

The term on the left hand side of is connected with the midpoint rule of the approx- 
imate integration, while the term on the right hand side is connected with the trapezoidal 
rule. Then the inequality (O can be regarded as an example of an extremality for quadra- 
ture operators. Many extremalities are known in the numerical analysis (cf. J4), cf. also 
and the references therein). The numerical analysts prove them using the suitable differen- 
tiability assumptions. As we will show in this paper, for convex functions of higher order 
some extremalities can be obtained without assumptions of this kind, using only the higher 



Date: February 2008. 

2000 Mathematics Subject Classification. Primary: 26D15, 41A55, 41A80; Secondary: 26A51, 65D30, 
65D32. 

Key words and phrases. Approximate integration, error bounds, extremalities, Hadamard-type inequalities, 
higher order convexity, quadrature rules, support theorems. 



2 



SZYMON WASOWICZ 



order convexity itself. The suppor-type properties play here the crucial role. A general the- 
orem of this nature was recently proved by the author in fill . The obtained extremalities 
are useful in proving error bounds of quadrature operators under regularity assumptions 
weaker from the commonly used. The results of this sort are also known, however our 
method seems to be quite easy. But the price we must pay is high: the obtained error 
bounds are far to be optimal (cf. 0). Some results concerning the inequalities between 
the quadrature operators and error bounds of quadrature rules, which are partial cases of 
our results, can be found in author's earlier papers lfT7l[l6l[l9l . The paper 1151 contains the 
extension of our results concerning convex functions to the functions of several variables. 

For n 6 N denote by n„ the space of all polynomials of degree at most n. Recall that 
a linear functional T defined on a linear space X of (not necessarily all) functions mapping 
some nonempty set into K is called positive if 

/ < g => T(f) < T(g) 

for any /, g € X. An important class of positive linear operators form the conical com- 
binations of the involved function at appropriately chosen points of a domain. Obviously, 
if a domain is a real interval, then quadrature operators with nonnegative coefficients are 
linear and positive. 

Dealing with a problem of approximate computation of the integral over an interval 
[a, b] it is enough to change the variable and to compute it over a fixed interval. The 
interval [—1,1] is frequently used. For a Riemann integrable function / : [— 1, 1] — * K let 

1(f) := j f(x)dx. 

Definition 2. Let T be a linear functional defined on a linear space of (not necessarily 
all) Riemann integrable functions mapping [—1,1] into M containing II n . We say that T is 

exact on II n if T(p) — I(p) for all p € H n . 



2. Convex functions of higher order 
Recall that the divided differences are defined as follows: [xi;/] := f(x±) and for 

k e N 

[x 2> ...,x k+1 ;f] - [x!,...,x k ;f] 
[xi,...,x k+ i;f\ ■■= . 

Xfc+1 - Xl 

Definition 3. Let / C K be an interval and n € N. A function / : / — * K is n-convex if 
[xi, . . . , x n +2] f] > for any distinct xi, . . . , x n +2 € /. 

There is an easy to imagine geometrical equivalent condition of n-convexity (for the 
proof cf. e.g. 0Q31)- 

Proposition 4. A function f is n-convex if and only if for any n + 1 distinct points 
xi, . . . , x n +i G /, the graph of the (unique) polynomial p € Il„ interpolating f at these 
points, passing through each point (xi, /(xj)J, i = 1, . . . , n + 1, changes the side of the 
graph of '/ (always p ^ f on [x n ,X n +i]). 

Then trivially 1 -convexity reduces to the classical convexity. 

Convex functions of higher order are very well known and investigated (see e.g. 
[lOl [TT1 FPU ). But up to now there is no common terminology, which sometimes may be 
confusing. 
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The first person who dealt with the topic in question was Hopf. He considered in his 
dissertation [ 7 1 from 1 926 the functions with nonnegative divided differences without nam- 
ing them at all. The notion of higher order convexity was introduced by Popoviciu in his 
famous dissertation [11 1 from 1934 exactly in the sense of the above Definition [3] The 
Kuczma's monograph [9| devoted to functional equations and inequalities in several vari- 
ables as well as the classical Roberts and Varberg's book on convex functions fl"3l use the 
same terminology (according to which an ordinary convex function is 1 -convex). During 
the years another way of naming convex functions of higher order became popular. Some 
authors (cf. e.g. [TOj) call a function / to be n-convex if \x\ , . . . , x n +i ; /] ^ (then 
a convex function is 2-convex). Now these two terminologies appear simultaneously in 
the literature. The first one is concentrated on the maximal degree of the interpolating 
polynomial, while the accent of the second one is put on the dimension of the space of 
polynomials of degree not exceeding some natural number. Then the second terminology 
is more coherent with convexity with respect to Chebyshev systems (cf. [8 1). Both conven- 
tions have some advantages and disadvantages and it is not the author's intention to judge 
neither which one is better nor which one is classical. 

Having in mind the above remarks let us declare that in this paper we understand the 
higher order convexity in the sense of Definition [3] 

Convex functions of higher order have many regularity properties. For details see j9] 
[TT][T3]. The paper ifTUl contains a brief survey of the topic given in one place. Below we 
list the properties which either we use in the paper or we discuss below. 

Theorem 5. If f : [a, b] — > M. is n-convex then f is continuous on (a, b) and bounded on 
[a, b}. 

Corollary 6. Iff : [a, b] — ► M is n-convex then f is Riemann integrable. 

Theorem 7. The real function f defined on an open interval I is n-convex if and only if 
is convex on I. 

Corollary 8. If the real function f defined on an open interval I is n-convex then 
/i_ exist on I and f^ exists almost everywhere on I. 

Notice that there are n-convex functions which are not n times differentiable (e.g. 
f(x) — \x\ for n = 1). Sometimes what is proved under differentiability assumptions, 
holds in fact for any n-convex function without further assumptions. We return to this mat- 
ter in Section [3] However, the following result requiring the differentiability assumption 
seems to be important (cf. j9j[TTl[T3|, for a quick reference cf. also [16 Theorems A and 
B] and Ql Theorem D]). 

Theorem 9. Assume that / : [a, b] — > M is (n + \)-times differentiable on (a, b) and 
continuous on [a, b]. Then f is n-convex if and only if (x) ^ 0, x € (a, b). 

It is well known that a convex function defined on a real interval admits an affine support 
at every interior point of a domain. In the paper lfl8l we have proved a general suppor-type 
result for convex functions of higher order. Four special cases ( ifTHl Corollaries 8-11]) play 
the crucial role in the proofs presented in this paper. 

Theorem 10. /// : [a, b] — * R is (2n — \)-convex and x\ 1 . . . , x n € (a, b), then there 
exists a polynomial p € n2n_i such that p{xi) — f(xi), i = 1, . . . , n, and p ^ f on [a, b\. 

Theorem 11. If f : [a, b] — > K is (2n — l)-convex and x\ = a, X2, ■ ■ ■ ,x n S (a, b), 
x n+ \ = b, then there exists a polynomial p € n2„_i such that p(xi) = f(xi), i — 
1, . . . , n + 1, and p ^ f on [a, b}. 
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Theorem 12. If f : [a, b] — * R is In-convex, x\ = a, X2, ■ ■ • , aVi+i G ( a i frten f«ere 
exists a polynomial p £ II 2n rac/i that p{xi) — f(xi), i = 1, . . . , n + 1, ant/ p ^ f on 
[a,b]. 

Theorem 13. Tjf / : [a, b] — > R is 2n-convex, a;x, . . . ,x n £ (a, 6) an<i x n+ i = b, then 
there exists a polynomial p £ n 2 „ such that p{xi) = f(xj), i = 1, . . . , n + 1, and p ^ / 
on [a, b]. 



3. EXTREMALITIES FOR QUADRATURE OPERATORS 

What we recall below is very well known from the numerical analysis (cf. e.g. JT] 
|6] [121 [141 [20] EU 1221). Let P n be the n-th degree member of the sequence of Legendre 
polynomials. 

Gaus-Legendre quadratures. For / : [— 1, 1] — ► R and n £ N let 

n 

Gn(f) ■= ^2mf(xi), 

1=1 

where x\, . . . , x n are the roots of P n (which are real, distinct and belong to (—1, 1)) and 



2(1-*?) 
(n + l) 2 P„ 2 +1 (x i ) : 



W i — i 1 1 \ 2 p2 t — \ ' * — 1, . . . , 71. 



Then £ n is exact on n 2 „_i. If / G C 2 "([-l, 1]) then 

(3) ^=^ + (2n^)[W /(aW)( ° 
for some £ G (— 1, 1). 

Lobatto quadratures. For / : [-1, 1] -> R let £ 2 (/) := /(-I) + /(l). For n G N, 
n ^ 3, let 

n-1 

£„(/) := + w n f(l) + Wif(xi), 

i=2 

where x 2 , . . . ,x n -i are the roots of P^-i (which are also real, distinct and belong to 
(-1,1)) and 

2 2 

W\ = W n = — r, Wi — — r t ; -, % = 2, .... 71 — 1. 

n(n-l)' 7i(7i-l)p2_ l(a;i )' 
Then £„ is exact on U 2n ^ 3 . If f £ C 2 "" 2 ([-l, 1]) then 

(4) 1(f)- C (f)- "("-D 3 *- 1 ^)!] 4 f (2n - 2) (0 
(4) ±U)-L n U) ( 27l _i)[(2n-2)!] 3 ; Uj 

for some £ G (—1, 1). 
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Radau quadratures. For / : [— 1, 1] — ► R and n e N, n > 2, let 

<(/) :=^/(-l) + X]^/(^): 

i=2 

where x%,...,x n are the roots of the polynomial 

P n _ 1 (a;)+P n (a;) 

Qn-l[x) = — 

X + 1 

(again real, distinct and belonging to (—1, 1)) and 
2 1 

»1 = o ) Wi 



(i-^)K-i(^)] 2 ' 

Then 7^ is exact on II 2 „_ 2 . If / e C 2n_1 ([-1, 1]) then 

for some £ £ (—1, 1). 

In (18 1 we also considered the operator 

K(/) :=<(/(-•))■ 

It was, in fact, defined in terms of orthogonal polynomials. However, these two definitions 
coincide. This is not difficult to check. We would not like to go into details since this 
is not the goal of the paper. Let us only mention that (changing the way of naming and 
numbering the abscissas and weights) we have 

n 

and TV n is exact on Tl2n-2- The error term of lU n is similar to ©, precisely, if / £ 
C 2n_1 ([-l,l])flien 

(6) J (/) = ^(/)- 22n [ (2ni%3 )!]4 / (2 "- 1) W 

for some 77 £ (—1, 1). 

As we can see, all the weights of the above quadratures are positive, so these operators 
are positive. This is also the case for many other quadratures. However, there are the 
quadratures with negative coefficients (e.g. among the Newton-Cotes formulas). 

Now we can prove the main results of this section. 

Theorem 14. Fix n £ N. Let T be the positive linear operator defined (at least) on 
a linear subspace o/R[ _1,1 l generated by a cone of (2n — V)-convex functions. Assume 
that T is exact on \\2n-i- If a function f : [— 1, 1] — ► R is (2n — l)-convex then 

(7) Q n (J) < T(/) < £„ + i(/). 

Proof. By Theorems [TOl and [TTI there exist two polynomials p,q £ H2n-i interpolating / 
at the abscissas of the operators Q n , C n +i, respectively, such that p ^ / ^ q on [—1, 1]. 
Since Q n = £ n+1 = I on n 2n -i, we get 

Qn{f) = Gn(p) - T(p) = T(p) < T{f) < T(q) = l(q) = £ n+1 (q) = £„ +1 (/). □ 
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Theorem 15. Fix n G N. Let T be the positive linear operator defined (at least) on 
a linear subspace o/Rl -1,1 ] generated by a cone of In- convex functions. Assume that T 
is exact on IT^n- If a function f : [— 1, 1] — > R is In-convex then 

(8) ti(/)a(/)a; +1 (/). 

Proof. Use Theorems [T2l and IT~3l and the abscissas of the operators 7Z l n+1 , TZ r n+1 , respec- 
tively, and argue similarly as in the proof of Theorem[l4] □ 

We would like to emphasize two particular cases of the above results. The first one 
concerns the inequalities of Hadamard type. The assertions of the Corollary below were 
proved in |[T8l Propositions 12 and 13] (cf. also the earlier paper [1|, where these results 
were obtained by another method). 

Corollary 16. Fix n E K If f : [— 1, 1] — > R is (2n - \)-convex then Q n (f) < 1(f) < 
Iff ■ [-1, 1] -» R is In-convex then K l n+1 (f) < 1(f) < 

Proo/ Use Theorems [H and [B] for T=l. □ 

The second important case is connected with quadrature operators. 

Corollary 17. Fix n, N £ N, £i, . . . ,£n G [-1, 1] and Ai, . . . , Ajv ^ 0. Lef 

/or/:[-l,l]^R. 

i=l 

(i) Tjf T exacf on Il2 n -i, fnen Q n (f) ^ ^"(/) ^ <Cn+i(/) /or any (2n — \)-convex 
function f : [— 1, 1] — > R. 

(ii) TfT ;s exact on IT^n, then lZ l n+1 (f) $J ^~(/) ^ T^-n+iif) f or an y convex function 
/:[-l,l]-K. 

Proof. The operator T trivially fulfils the assumptions of Theorems [14] and [15] respec- 
tively. □ 

In the numerical analysis the inequalities the above type are called extremalities. The 
extremalities of Corollary [17] (Qi were earlier proved in [4 Theorem 6] under the assump- 
tion of 2n-times differentiability. The proof given there is based on taking double nodes. 
The author independently used in |[T8l exactly the same idea to prove support-type results 
of Corollaries 8-11 (quoted here in Theorems [TOlfOl l with no use of any differentiability 
assumptions. Thus, as we can see from the proof of Theorems [T4l and [TBI the extremalities 
in question are proved with no further assumptions, except higher order convexity itself. 
So, our results are more general than these of (4). 

We underline that the inequalities of Corollary [T7] do not hold for any quadrature oper- 
ator T. The exactness assumption (i.e. T = X for polynomials of the appropriate degree) 
is essential. 

Example 18. We have 

«/> = /(-f)+/(f 



<«(/) = |/(o) + 1 



5 / V 5 



n 3 (f) = -/(-i) + -^-/ f — b — 1 + -^g-/ 
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(cf. e.g. 1201 221 ). whence ^(exp) > £/2(exp) and "^(exp) > ^(exp). The exponential 
function is convex of any order (cf. Theorem |9). Let T = Q2- Then the inequality of 
Corollary [17] (0) does not hold for n = 3, and that of ([n} is not true for n = 2. Notice that 
for p(x) — x A we have Q2 (p) ^ I{p), so G2 I both on 1T5 and on II 



4- 



Now the question arises if there are other, i.e. non-quadrature, operators approximat- 
ing the integral, for which Theorems [14] and Q3] are applicable. The positive answer given 
below shows that the extremalities for quadrature operators are special cases of more gen- 
eral inequalities. Not the form of the operators considered (linear combination, integral 
and so on) is important but two things play the key role: positiveness and exactness for 
polynomials of the appropriate degree. 

Example 19. For a Riemann-integrable function / : [—1,1] — > K let 

T(f) [/(-I) + /(1)] + TI J ,/(*)*• 

Then T is a positive linear operator exact on II3. Using Theorem [14] for n = 2 we obtain 
&(/) < T (/) < £ 3(/) for a 3-convex function / : [-1, 1] -> M. 

For the other non-quadrature operators approximating the integral cf. e.g. Q. 
4. Error bounds of quadrature operators 

In this section we show that the extremalities of Corollary [l7]may be applied to obtain 
the error bounds of the involved quadrature operator T using the regularity assumptions 
weaker from the commonly used. The results of this type are known in the numerical 
analysis. We would like to point that the inequalities of Hadamard type may be used in 
the approximate integration. But the results obtained by our method deliver error bounds 
which are far to be optimal. This is a price we have to pay for simplicity. Error bounds 
obtained in the numerical analysis under assumptions used by us are much better (see J5] 
and the references therein). 

For /GCQ-1,1]) denote H/lloc := sup {\f(x)\ : x £ [-1,1]}. 

Lemma 20. Fix k G N, k ^ 2. Let IC, 7~ be linear operators defined on a linear subspace 
of R[ _1,1 l containing all the functions involved below with the following properties: 

(i) there exists an a > such thatl(f) < K(f) + a\\f^ \\oofor all f £ C fc ([-1, 1]); 

(ii) Iff : [-1, 1] — > M is (k - \)-convex then £(/) < T(f); 

(iii) T(p) —I{p) for p{x) — x k . 

Then \l(f) - T(/)| < 2a\\f^\\ 00 for any f £ C fe ([-1, 1]). 
Proof. By (01 and (Ej} we get 

(9) 2:(/)-n/)<«ll/ (fc) lloo 

for any (k — l)-convex function / £ C fc ([-1, 1]). 

For an arbitrary function / £ C k {[— 1, 1]) define now g(x) := ^ x k . Then = 
||/ (fe) ||oo, whence 1/(^1 < on [-1, 1], which implies {g-f) {k) >0and(p+/)W ^0 
on [—1, 1]. Therefore by Theorem[9]the functions g — f, g+f are (k— l)-convexon [—1, 1]. 
By the triangle inequality 

||(.9-/) (fe) ||c O ^2||/( fc )|| 0C and || (g + || ^ < 2||/W|| 0Q . 

Now we apply © to g — / and g + f. Then the desired inequality follows by linearity, the 
assumption (liTit and the above inequalities. □ 
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For n£N let 



4 n+1 (n!) 4 4" +1 (n + l)(n!) 



M4 



a 2n '■— -77: ; , u/n rPTo ; Oi2n+l 



(2n + l)[(2n)!] 3 ' ' [(2n+l)!] 3 

Theorem 21. F/x A; £ N, k 2. Let T be a positive linear operator defined on a 
domain as in Lemma l20l If T is exact on life, then \l(f) — T{f)\ ^ C-kWf^Woa for 
anyf£C k {[-l,l}). 

Proof. If k is even, k — 2n, then use LemmaEOlfor a = and JC — Q n . The condition 
(i) is fulfilled by (O, (m) holds by TheoremfBland (Hi) by the assumption. 

Similarly, if k is odd, fc = 2n + 1, then use Lemma l20l for a = " 2 ^ +1 and /C = 
^n+i- The condition (i) is fulfilled by (O, (ii) holds by Theorem [T5l and (in) by the 
assumption. □ 

In the assertion of Theorem [14] there is an operator £„+i on the right hand side of the 
inequality dT) and in the statement of Theorem Q3] there is an operator IZn+i on the right 
hand side of the inequality ([8). We could prove the result similar to Theoreml2"Tl involving 
(in the proof) these operators. However, the error bound obtained in this way will not 
improve that of Theorem[2TJ Namely, the absolute value of the constant of <j4j (take n + 1 
instead of n) is greater from the similar constant of For the operators Tl l n+l and lZ 1 n+1 
the absolute values of both constants of ©, (|6]l are the same. 

For the quadrature operators we immediately derive from Theorem|2T|the following 

Corollary 22. Fix k,N £ N, k > 2, £1, . . . ,£ N £ [-1, 1] and Xi, . . . , X N > 0. Let 

N 

T(f) 1=5^/(6) for f : [—1, 1] — > M. 

i=l 

If T is exact on U k , then \l(f)-T(f)\ < a fe 00 /or any / € C fe ([-l,l]). 

Using this result we will now give the error bounds of Gauss-Legendre, Lobatto and 
Radau quadratures under regularity assumptions weaker from the commonly used. Denote 
by |_-J the floor function, i.e. [^J = maxjfc £ Z : k ^ x}, 

Proposition 23. Let k, N £ N, k ^ 2 and N > [_§J. Iff £ C k ([-1, 1]) then \l(f) - 
G N (f)\ SjafcH/WlU 

Proo/ If TV > L|J then 2N - 1 ^ k, whence Q N is exact on II fc (cf. ®). Then the result 
follows immediately by Corollary [22] □ 

Proposition 24. Let k, N e N, k ^ 2 ant/ iV > [|J + 1- #7 £ C*([-l,l]) then 
\l(f)-£ N (f)\ ^a fe ||/ (fc) ||oo. 

Proof. This is also an immediate consequence of Corollary l22lsince JV > L|J + 1 implies 
2N - 3 > fc and then £at is exact on II fc (cf. ©). □ 

Proposition 25. Let k, N £ N, k ^ 2 one/ TV > [^J- /// € C fc ([-1, 1]) f/ze« |^(/) ~ 
^5v(/)| ^ ^fcll/^^lloo- The same assertion holds for the operator lZ r N . 



Proof. By TV > [^J we get 27V - 2 > ft and we can see (cf. ©, ©) that ft^ and ft 
are exact on IT^, which, together with Corollary[22j concludes the proof. □ 



JV 
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